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Abstract 

This paper describes an approach for controlling the at- 
titude of statically unsable thrust-levitated vehicles in 
hover or slow translation. The large thrust vector that 
characterizes such vehicles can be modulated to provide 
control forces and moments to the airframe, but such 
modulation is accompanied by significant unsteady flow 
effects. These effects are difficult to model, and can com- 
promise the practical value of thrust vectoring in closed- 
loop attitude stability, even if the thrust vectoring ma- 
chinery has sufficient bandwidth for stabilization. The 
stabilization approach described in this paper is based on 
using internal angular momentum transfer devices for sta- 
bility, augmented by thrust vectoring for trim and other 
“outer loop” control functions. The three main compo- 
nents of this approach are: (1) a z-body axis angular mo- 
mentum bias enhances static attitude stability, reducing 
the amount of control activity needed for stabilization, 
(2) optionally, gimbaled reaction wheels provide high- 
bandwidth control torques for additional stabilization, or 
agility, and (3) the resulting strongly coupled system dy- 
namics are controlled by a multivariable controller. A 
flight test vehicle is described, and nonlinear simulation 
results are provided that demonstrate the efficacy of the 
approach. 

1 Introduction 

1.1 Background 

There has been recent interest in flight vehicles, capable 
of VTOL and hover, that are statically unstable or neu- 
trally stable [1]. Examples of such vehicles include small 
observation platforms [2], [3] in which an optical sensor 
is mounted above a levitation-producing ducted fan, and 
rescue vehicles [4] consisting of a payload-bearing plat- 
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form with a number of ducted fans distributed symmetri- 
cally about its perimeter in the plane of the platform. In 
the state of the art, this type of vehicle is controlled by 
thrust vector variation, possibly with additional aerody- 
namic features implemented to enhance stability in the 
presence of gusts [2], [3]. 

Using thrust vectoring or differential thrusting to sta- 
bilize a statically unstable airframe poses two related dif- 
ficulties. First, because of the airframe’s lack of stabil- 
ity, the controller must operate with high authority and 
bandwidth in order to contain and suppress excursions 
in attitude. Accurate system models are generally re- 
quired to successfully design such controllers. The second 
difficulty is that such models are not available. This is 
because the active thrust vector variations required for 
stabilization in the presence of disturbances induces un- 
steady flow phenomena that distort the control command 
in a complicated manner that currently defies compre- 
hensive modelling (see for example [5]. These difficul- 
ties come together when the control engineer attempts to 
design the required high-performance controller while ig- 
noring or under-modelling the effect of the unsteady flow 
phenomena on the closed-loop dynamics These difficul- 
ties become more evident for smaller inertia VTOL vehi- 
cles because they respond more rapidly to disturbances 
so that correspondingly more rapid feedback control is 
necessary for vehicle stabilization. 

1.2 A New Approach 

This paper describes the configuration, dynamics, and 
control of a statically unstable flight vehicle that levi- 
tates by means of thrusters with thrust vectoring devices, 
but augments the thrust vector control with internal an- 
gular momentum transfers to enhance attitude stability 
and provide secure command tracking performance in the 
presence of external disturbances, such as wind gusts and 
in-flight payload variations. 

The salient features of such a vehicle are: 

• The vertical component of the thrust vector levitates 
the vehicle. 

• The thrust vector is produced by a plurality of thrust 
generators, such as ducted fans. The direction of the 
thrust vector is modified by changing the orientation 


of the thrust generators, such as by gimballing, or by 
diverting their flow. The contributions from the plu- 
rality of thrusters produce net forces and moments 
about the airframe’s center of mass that can be used 
to trim the vehicle and for control. For the exam- 
ple system explored further in the sequel, thrust is 
provided by four ducted fans, each with axis of rota- 
tion fixed parallel to the z body axis, and the thrust 
direction is modulated by a pair of thrust vectoring 
vanes under each fan. 

• Static attitude stability is enhanced by an angular 
momentum bias in the z body axis. This bias is pro- 
duced by a plurality of rotating masses, such as ded- 
icated momentum wheels. Alternatively, when em- 
ploying fans or impellers whose rotation is aligned 
with the z body axis for levitation thrust, the z- 
body angular momentum bias can be generated by 
augmenting their moments of inertia about their ro- 
tation axes. This is done by redistributing mass to- 
ward, or introducing additional mass near the fan or 
impeller radius for one or more of the fans or im- 
pellers. An example of this is attaching a ring to the 
tips of a levitation fan. 

• Additional attitude control authority is provided by 
a plurality of control moment gyros (CMGs). These 
interact with the thrust vector control in jointly 
generating moments for attitude stabilization; the 
thrust vectoring system is also used to desaturate 
the CMGs. 

• The control effectors described above are jointly 
commanded by a control system for stabilization 
and command tracking, and for desaturation of the 
CMGs. 

The advantages of the above vehicle are: 

• The vehicle’s open-loop attitude stability is enhanced 
by the z-body axis angular momentum bias. This 
reduces the authority required from the active stabi- 
lization control system, and permits a reduction in 
the bandwidth of the thrust vectoring control com- 
mands. 

• It is straightforward to develop accurate control de- 
sign models of the devices for internal angular mo- 
mentum exchanges. Because of this, they can be 
used in high-bandwidth, high-authority control sys- 
tems. This permits enhanced disturbance rejections 
and maneuverability. 

• The use of internal angular momentum exchange de- 
vices relaxes the performance requirements on the 
thrust vector control effectors, since their role can be 
reduced to as little as that of providing trim, com- 
mand tracking, and CMG desaturation. 

Figure 1 illustrates the basic differences between the 
state-of-the-art approach and our new approach. 
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Figure 1: Inner-loop attitude control by differential and 
thrust vectoring, state-of-the-art (upper), and by internal 
torques aided by bias momentum, our approach (lower). 


1.3 Outline of Paper 

This paper outlines recent results in our continuing R&D 
efforts as first reported in [6]. In the next section, a re- 
search vehicle under development at NASA Langley Re- 
search Center for maturing this control technology is de- 
scribed, and its equations of motion are given. Section 3 
discuss its flight dynamics and control aspects including 
vehicle trim, gyric stability, and control structure. The 
sizing of the bias momentum for this new class of “hover- 
ing spacecraft” vehicles is also discussed. In section 4, we 
describe a current flying test platform system under de- 
velopment and outline its simulation-based performance 
assessment and comparison with a similar vehicle not us- 
ing momentum exchange devices. Section 5 provides con- 
clusions. 


2 Equations of Motion 

To derive the equations of motion for a generic NFTP 
system 1 , a Newton-Euler formulation is used, and in par- 
ticular, the sequence of formulation and notation found 
in [7] , and extended to multiple gimbaled spinning bodies 
with vane rotational motion for vector thrusting. Only a 
summary of the results are given here and the interested 
reader is referred to [8] for details. 

Figure 2 shows a schematic of the NFTP describing 
its main components. 


1 “Generic” in the sense of a flying platform propelled by multiple 
ducted-fan and vane system and augmented with multiple wheel 
system 
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Figure 2: Schematic of NFTP. 


R Rigid platform 
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Ob origin of platform fixed axes 
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0 Cb c/m of fl 

Of, origin of i-th ducted fan subsystem axes 
0 Gi origin of i-th CMG frame 

Let uj Si a,i, f Ib&, and r/i/, denote the relative angular 
velocities of the zth fan, momentum wheel, and the gimbal 
system relative to platform, respectively. Additionally, 
let iliSi), denote the relative angular velocity of the zth 
CMG wheel relative to its gimbal frame. For the CMG, 
assume that the center of mass of rotor is fixed in F b , 
and that the gimbal axis, g , is fixed in Fb, while the spin 
and output axes are defined as s and d respectively. For 
more details on CMGs, see for example [9]. We define the 
following frames: 
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In summary, we consider an 19 degrees of freedom NFTP 
model which includes the following: 

3 positions of Ob relative to O 

3 rotations of R relative to F 0 

4 rotations of fan relative to Fb 
1 rotation of B relative to Fb 

4 rotations of CMG gimbals relative to Fb 
4 rotations of CMG wheels relative to F G . 


The following momenta variables are used in describ- 
ing the equations of motion: 


p = m_u + ii/ x (1) 

4 

x I^ F .LO Si ai 

i - 1 
4 

+ y + IwfliSi) + Ib^bZ (2) 

i-1 

h ai = dp. zb • a; + I^p.u! Si , i = 1, ... ,4 (3) 
h Bz = I%z-(u+n B z) (4) 

h Wil = F^« + fti), z = l,...,4 (5) 

hCi 2 = +Vi): i = 4 (6) 


The vectors, and denote the platform translational 
and rotational inertial velocities, respectively, while 
denotes the first moment of inertia about Ob for the sys- 
tem having a mass to. The vectors, p and h^, denote the 
the translational and rotational momentum of the sys- 
tem, respectively. The angular velocity term, ujJ‘ where 
i = 1, . . . , 4 and j = 1,2,3 (corresponding to Sj, gi, and 
di), denotes the platform angular velocity components in 
the jth axis of the zthe gimbal frame. In particular, if we 
represent the platform angular velocity in terms of body 
frame so that = F^u>, the above angular velocity terms 
can be expressed in matrix form as 


ujJ = Kui, 0 J 2 = Ff 0 a;, W 3 = Gu> (7) 

where 
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and for z = 1, . . . , 4 


Kvi) - 

[ cos r]i 

0 

- sin T)i ] Cq. 

( 8 ) 

k Qi = 

[0 1 

0 ] 

(~10 
0 Gi 

(9) 

aim) - 

[ sin r]i 

0 

cos Pi ] Cq. 

( 10 ) 


For practical reasons, the rotational momentum is 
defined about the point Ob, which is the center of the 
platform frame, Fb, but need not be the center of mass 
of S due to practical limitations and/or not accurately 
known variable mass payload and its distribution during 
flight. Naturally additional complications arise due to the 
fact that reference point Ob is neither the center of mass 
nor fixed to an inertial frame. The remaining variables, 
h ai ,h Bz , , hc i2 , are absolute rotational momenta vec- 
tor components of the zth rotating fan, momentum wheel, 
CMG wheel, and CMG system, along their respective pri- 
mary spin axes (i.e. a,i,z,§i,gi). As reference points for 
the above momenta components, their respective center 
of masses are used. Of course the center of masses for 
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these actuator subsystem hardware can be accurately de- 
termined and are not expected to vary during operation. 

With the above choice of variables and notations, the 
equations of motion are summarized as follows (for details 


[8]): 

4 


P 

= mg + V T. + / +/ 

(11) 

—• > 

^ „ 1 *jet *aero 

i=l J 
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A 

= — jj>xp+_c^xg + r -(- 

i—1 

x T ] 


+ T + T 

*jet >aero 

(12) 

h ai 

= T Fi 

(13) 

h'Bz 

= T B z 

(14) 

hw ix 

= Tsi 

(15) 

h Ci2 

= [(iCj — ^Ci)bJ l‘ + Iw^i] + Tgi 

(16) 


where i = 1, ... ,4, i.e., 4 sets of ducted fans and sin- 
gle gimbal variable speed CMGs are assumed, but clearly 
any number can be used in these equations. The terms, 
T = T +T , t = t + t , denote external aero- 

— H — >Wi ~*Wi — >b; 

propulsive net thrust and torques on the system due to 
vector thrusting and differential throttling of the ducted 
fan and vane subsystem. The terms / , / , r , 

— >jet — >aero J e ^ 

and t denotes the aerodynamic loads on the system 

— >aero 

due to significant vehicle airspeed. In this paper, we focus 
on the subclass of problems where the vehicle is flying at 
low airspeeds, due to its particular difficulty in attitude 
stabilization. Hence we do not include these terms in this 
paper but is addressed in [8]. The terms t§. and Tg. de- 
notes the internal torques on the i-th CMG wheel due to 
the wheel motor, the internal torques on the i-th CMG 
subsystem due to the gimbal motor, respectively. The 
variable t b denotes the internal torque on the momen- 
tum wheel due to its motor, while t f denotes net torque 
on the i-th fan due to a combination of external fan rota- 
tional drag and internal fan motor torque, including back 
EMF effects. 


2.1 Aerodynamic and propulsive loads 

Figure 3 shows a set of frames used to conveniently de- 
scribe various loads produced by engine thrust, vane de- 
flections, and fan rotational drag. 


2.1.1 Net propulsive forces acting on Op i 

T.(u Si ,0i) = F^C^T^Oi), Ti G U 3xl (17) 
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Figure 3: Reference frames for modeling loads produced 
by Fans, Vanes, and CMG. 


2.1.2 Net propulsive moments acting about Op i 

_r .(u> Si ,9i) = F £ Cp.Ti{bJ Si ,9i), n G 77 3x1 (19) 

where 

l Xi \ 

hi 9iU 2 Si '■= fi{9i)u>si 

A J 

(20) 



2.1.3 Rotational Drag on i-th Fan along a; 

CH ' 1> W . = ~ T Wi Os,)> T Wi € ft ( 21 ) 

where 

T Wi (Usi) « C R^H ( 22 ) 

Fan rotational drag only in a; axis (i.e. excludes vane 
deflection effects). 


2.1.4 Torques on i-th Fan motor along a; 


Oi ’ = T ai (W fl 


,0J 


cmd 


), 


G 77, i = 1, ... ,4 (23) 


where 

Ta; (w s< , LO c ™ d ) ss k EMF . uj Si + k 3 u >^ d (24) 

The torque modeled in equation 24 is a combination of 
back EMF and armature current torques. The net torque 
on the fan system along spin axis due to rotational drag 
and motor 
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2.2 Scalar form 


For modeling and computational purposes it is more con- 
venient to express the motion equations in scalar/matrix 
form. In particular, we choose to express vector quan- 
tities in terms of their components in F b , Fp . , and Fq : 
frames: 
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F Iv 
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where all column matrix variables belong to 1Z 3 x 1 except 
J G 7Z 3x3 . These quantities include geometrical and iner- 
tial parameters of the vehicle that are constants in body 
frame, and state variables and forces and moments whose 
components can vary in body, gimbal, or fan frames. 

The first order momenta equations for the total sys- 
tem are summarized as follows 
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P 
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0 

X 

3 

1 

/ p X , 

ml 3x3 

— V X — (jJ x 

1 >> j + 

c x 


C t f 0 1 f T(u s ,9) \ 

C bF C t f \ \ t(uj s ,0) J 


(26) 


through the effects of gravity. Its inadvertent effects on 
rotation is due to the limited knowledge in the exact 
location of the system center of gravity. If the center 
of gravity of the whole system is known, then, O b can 
be chosen as such and the rotational and translational 
motions uncouple. 

In the above equations, the superscript (-) x denotes 
the skew-symmetric matrix associated with vector cross 


products, for example, u x = 


0 — u 3 LO2 

U 3 0 —LUi 

—L02 u>i 0 


For the definition of augmented states and variables in- 
cluding uj s ,h a ,D,r],h Wl ,hc . 2 > -T T , T w , T a , Tg , Tg , T F , 

and augmented parameters including 

C F , C bF , A, , If, If, J“ F , If,, If, see [ 8 ] for details. 

For kinematics, we denote the position of the ori- 
gin, O b , of the platform-fixed frame, F b , with respect 

to inertial frame, F 0 as R = F so that the iner- 

y( - >b 

tial velocity of O b is = Ff'v = and the iner- 

tial velocity components in F a frame can be written as 
£ = C((3) T v where we choose to parameterize the direc- 
tion cosine matrix relating F a to F b by Euler Parameters, 
/ 3 := (/ 3 o,/ 3 i,/ 32 ,/? 3 ) T € lZ ixl . For the form of this direc- 
tion cosine matrix parameterization and the time deriva- 
tive of these Euler parameters, see for example [10]. 


3 Flight Dynamics &; Control 


and for ducted fan along spin axis are 


3.1 Trim 


h a =Tp{L0 s ,uf nd ) (27) 

and for momentum wheel along spin axis, z 

h Bz =T Bm (n B ,QT d ) ( 28 ) 

and for CMGs along spin and gimbal axes are 

f h Wl 1 _ f Tg X 

l he. 2 J \ — diag(Go)) [(If — If)Kio + /“^O] + Tg J 

(29) 

The momenta is related to the velocities as follows 
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where M = 


04x3 

IfA T 

0lx3 

(0 0 If) 

04x3 
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00 
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ml 3x 3 


—c 
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, and T 4 diag (/“ F , If, If,, If). 


In the above choice of dynamical equation coordinates, 
the inertial attitude of the platform (i.e. C) enters 


We initially consider the flight dynamics as referenced 
to trim conditions which we define as attitude hold while 
moving at constant translational speeds. This of course is 
suited only for cruise and hovering flight segments. The 
trim conditions are found by calculating any combina- 
tions of four vane angles and four fan speeds necessary to 
hold a given desired attitude, at a constant translational 
velocity. With variable speed CMGs we choose a nom- 
inal wheel speeds and zero gimbal rotation rate as the 
trim value. In summary, the trim conditions are defined 
by setting all momenta state derivatives, given in equa- 
tions 26 to 29, to zero. In addition, we consider u = 03 X i 
to hold platform attitude to (3. 

Assuming the above conditions, the first order mo- 
menta equations for trim reduces to a coupled set of alge- 
braic equations [ 8 ]. We make the following observations 
from the trim conditions: 

• Because Cj = 0, p = mv, i.e. at trim, the linear 
momenta do not affect rotational equilibrium, i.e., 
v x p = 0 . 

• For a given steady or near trim flight condition, a 
combination of vector thrusting via fan speeds, uj s , 
and vane angle deflections, 9, are used to satisfy trim 
conditions. In the case when the vehicle is required 
to trim at airspeeds that result in significant vehicle 
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aerodynamic loads, as described in equations 11 and 
12, the above trim conditions can be slightly modi- 
fied. 

• For the ducted fan trim h a = 0, the fan control mo- 
ments, f a , exactly balance the fan rotational drag, 
t w as given in equation (25). The trim command 
can be computed once the required fan speed ui s is 
known from system momenta equilibrium equations. 

• For the momentum wheel trim, hs z =0, the wheel 
control moments, T™ otor exactly balances the wheel 
friction, T f ncUon _ The desired level of bias momen- 
tum for the overall system defines the momentum 
wheel trim speed, Clg, which will determine its corre- 
sponding command, 

• For each CMG wheel trim where hw.t = 0, the CMG 
wheel motor must exactly balance the CMG wheel 
friction at a chosen trim speed Q. 

• For each CMG gimbal system trim where hc. 2 = 
the CMG gimbal motor torque must be zero corre- 
sponding to gimbal rate trim, rj = 0. 


The trim conditions can be expressed only in terms of the 
trim angles, 9 , and fan angular velocities, u> s , by using the 
vector thrusting models as given in equations (17) to (24) 
for the forces and moments produced by fan rotation and 
vane deflections. These key subset of trim conditions can 
be rearraged to the following compact form: 

+ (j>0) = 0 6xl (31) 


where 

M{9) 

m 


\ c? 

[ h C »F 


ml 3x3 

l „x 


E(6) + 

C0)g o 



F(9) (32) 


(33) 


define the geometry, kinematics, mass properties, sec- 
ondary trim actuators, and vehicle aerodynamics. The 
scalar constant b 0 meters is a reference moment arm to 
help scale the combined force-moment equations to force 
units, which helps in the selection of a treshold numerical 
value for defining equilibrium. 

For any given set of variables, (/3, 9 , Cj s ,v), that satis- 
fies the vector thrusting trim equations (31), one can de- 
duce the corresponding ducted fan speed trim command 
ujf md . The remaining trim commands, momentum wheel 
speed Cl c ™ d , CMG wheel speed £l cmd , and CMG gim- 
bal rate ■q cmd , will clearly depend on their corresponding 
desired trim conditions (T2 ^ , T2 , 77 = 0), and their corre- 
sponding motor dynamics. Note that these details involve 
motor commands to regulate and track armature rotation 
speeds, and involve technology that is well understood 
and developed. 


Proposition 1 (Trimability) Given a platform atti- 
tude, (3, the system defined by equations (26)-(29) is 
trimmable, using vector thrusting ( using a combination of 
control vane deflections, 9, fan thrusts with fan rotation 
speeds Cj s , and reaction jets), if and only if there exists 9 
such that (j>0) £ Range[M(9)]. 

If the above trim conditions are met for (3 and 9, then, 
the corresponding set of all possible fan speed combina- 
tions can be parameterized as follows: 

cbg = — M (9) + <j>0) + if where if £ Null[M(0)] (34) 

Proposition 2 (Numerical Criteria for Trim) 

Given a platform attitude, f3, the system defined by 
equations (26)-(29) is trimmable, using vector thrusting 
if the following least squares equation error is satisfied: 

lk(/M)||2 :=\\[I-M(9)M(§)+]cf0)\\ 2 <e (35) 

The superscript symbol in (•)+ denotes its pseudo-inverse. 
Note that the treshold for zero equation error denoted as e 
should reflect practical limitations including (i) numerical 
error, in the order of machine epsilon, (ii) accuracy limi- 
tation in the actuation and measurement of vane angles, 
(iii)accuracy limitation in the fan thrust response to com- 
mand, including integer command resolution, and (iv) 
sensor accuracy limitaion in a high quality accelerometer. 
These factors are used to define an engineering treshold 
of “trim” for the NFTP during low airspeed operations. 
Since trim angles and speeds are non-unique if they ex- 
ists, we so choose the “smoothest” combination, see [8]. 

3.2 Gyric Stability about trim 

Figure 4 shows the perturbed dynamics about trim. In 



Figure 4: Perturbed dynamics about trim. 

general, the open loop dynamics involves complex per- 
turbed response due to gyroscopic coupling in the rota- 
tion axes and the coupling with translational degrees of 
freedom (in addition to coning due to axis misalignments 
and translational-rotational coupling due to off CG ef- 
fects). Hence for the sake of analytical tractability, we 
simplify by assuming a perfect knowledge of system CG 
location and the orientation of the system principal axes. 
With these assumptions, the tenth-order system charac- 
teristic equation for the linearly perturbed system shown 
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in Figure 4 can be written and solved explicitly. It turns 
out that the only non-zero eigenvalues have the form 


'VlO — ±/A 0 , 


where A 0 := 


\fJ1J2 


(36) 


Physically, this complex conjugate pair corresponds to 
the gyroscopic coupling in the pitch and roll axes. This 
oscillation or wobbling frequency given in equation 36 is 
analogous to the “precession frequency” in a dual-spin 
spacecraft (a gyrosat with a nonspinning carrier) as de- 
rived for example in Chapter 6 of [7]. Clearly, this wob- 
bling frequency increases linearly with the level of bias 
momentum and is inversely proportional to the rotational 
inertia of the vehicle. It can be shown that the transfer 
function matrix from external torques (zq, zq, 1/3) to small 
angular responses in roll, pitch, and yaw has the form 


5 p 1 
dp 2 
dps 


jG(s; hs, Ji, J2) 

02x1 

O 

X 

to 

1 

s 2 J 3 . 



where the roll-pitch transfer matrix is 


G(s; h B , Ji, J 2 ) — 



S 



JlJ2 


hs 

J1J2 

s 

J2 


(37) 


(38) 


Figure 5 shows the roll angle frequency response from dis- 
turbance torques in the roll and pitch axes. Notice that 
the linear yaw response depends only on the yaw distur- 
bance torque and is independent of roll and pitch motion. 
The main point is that bias momentum can significantly 
improve attitude stability by reducing the response to dis- 
turbances in roll and pitch. Furthermore, as the level of 
bias momentum increases, the roll response to roll distur- 
bances decreases more rapidly than to pitch disturbance 
torque. Alternately, the dynamically uncoupled system 
at zero bias momentum transforms to a system whose 
characteristic dynamics is gyrically dominated. 




Figure 5: Open loop frequency response in roll and pitch 
angles at different levels of bias momentum. 


3.3 Bias Momentum Sizing 

Consider the rotational dynamics of a generic dual-spin 
system with principal rotational inertias {J\,J 2 ^ 2 ) driven 
by disturbance torques (ti, 72, and T3) in all three axes. 
Suppose the main body rotates with angular velocity 


components in body coordinates denoted by (u>i,u> 2 , and 
W3) while the momentum wheel spins about the verti- 
cal (3-axis) at a fixed rate relative to the main body to 
produce a constant bias momentum oi h B . In addition, 
suppose the bias momentum in the wheel “dominates” in 
the following sense 

\h B \ » max (|( J3 - J2V3 1, |(J 3 - Ji)w 3 |) (39) 


Then, the dynamics of this generic dual spin system 
be approximated by 

can 

Jldll ~ 

-L0 2 h B + n 

(40) 

J 2 U 2 ~ 

wi h B + r 2 

(41) 

^3^3 = 

( J 2 - Ji)viU2 + t 3 

(42) 


These equations are not limited to “small” angular mo- 
tion (associated with linearizations) but rather to the 
dominance of the bias momentum, as defined in equation 
39. Although the yaw dynamics in equation 42 is nonlin- 
ear, the roll-pitch dynamics in equations 40 and 41 are 
linear. In fact this transfer function matrix for roll-pitch 
dynamics is identical (except 1/s term) to the dual-spin 
linearized roll-pitch dynamical equations 38. 

For this generic dual spin system, we consider a bias 
momentum sizing approach based on the angular velocity 
response of the main body to a set of disturbances. To 
obtain an analytical expression for roll-pitch response to 
random disturbance torques, we make some mild assump- 
tions on these unknown disturbances. Suppose the power 
spectral densities of the disturbance torques caused by 
unknown random turbulence, ri, and 72, with a common 
bandwidth of BW T , can be modeled as 

V (jv) = diag(p Tl , p T2 ), Mv < BW T (43) 


where 12 denotes the dummy frequency parameter. It can 
be shown that [8] the mean square response of the roll 
and pitch angular velocities can be expressed as follows 


£[w?+wf] 


y / 2max(p Tl ,p T2 ) 

V J\ J2\h B \ 


(44) 


where the above estimate applies for responses at lower 
frequencies v < min{BW T , A 0 ). Figure 6 shows the effect 
of scaling bias momentum on the worst case angular ve- 
locity response in the roll and pitch axis, in terms of the 
maximum singular value frequency response. It follows 
then that if we want the RMS angular velocity response 
to be bounded by a certain maximum desirable value, 
then the bias momentum can be sized to satisfy an in- 
equality based on equation 44. In addition to the above 
requirement on the RMS response of the angular veloc- 
ities, the sizing of the bias momentum should be such 
that the precession frequency lies outside the bandwidth 
of the disturbance, BW T . This bandwidth requirement 
can be satisfied if we choose bias momentum level to sat- 
isfy, based on equation 36, the following condition: 

A 0 := -j== > BW T (45) 

V ^1^2 


7 



Figure 6: Maximum singular value frequency response in 
roll and pitch anglular velocities for sizing bias momenta. 

3.4 Control system framework 

Figure 7 shows a flight dynamics and control schematic 
for a general NFTP vehicle with the use of both bias 
momentum wheel and a set of variable speed CMG’s. 
The vector thrusting commands include four fan speeds 



Figure 7: Flight dynamics and control schematic. 

(«J md ) and four vane angles ( Q cmd ) control platform at- 

titude (/ 3 ), and velocities ( v and uj). These vector thrust- 
ing commands for the fans and vane angles are primarily 
used for generating trim conditions and for open loop 
command tracking. For the general configuration with 
CMGs, a secondary but nevertheless crucial function of 
the vector thrusting commands will be to stabilize the 
CMGs so that their saturation tendencies will be miti- 
gated. To improve performance, fan speed regulators and 
vane servos are used. 

Notice that the Ducted Fan/ Vane aerodynamics and 
gravity directly influence system momenta p , h, while the 
CMG Dynamics has no influence on system momenta. 


Rather, the CMGs internally redistribute angular mo- 
mentum to control the platform angular velocity u and its 
attitude (3. A key advantage is that the CMG can gener- 
ate internal control torques to effectively and very reliably 
redistribute angular momentum to control the platform, 
independently of Ducted Fan/ Vane aerodynamics. This 
of course mitigates the uncertain but significant effects 
of unsteady aerodynamics induced by rapid control sur- 
face motions for attitude stabilization. This property of 
this control approach will particularly benefit smaller ve- 
hicles during hovering or operations at low airspeeds un- 
der significant wind turbulence because they will respond 
rapidly due to their smaller rotational inertias and will 
therefore require control effectors that can reliably gener- 
ate control forces and moments at higher bandwidth for 
attitude stabilization. 

While the bias momentum wheel is regulated to a cer- 
tain constant speed, fls, to provide directional stability 
in open loop, the CMG subsystem will consist of closed 
loop torque commands to gimbal and wheel motors for 
variable speed CMGs. These commands will be generated 
by a control law which will integrate the CMGs with the 
vector thrusting control system for a given momentum 
wheel augmented platform. In particular, an approach 
with extends the steering laws outlined in [16], [17], and 
[15] to account for factors including, trim requirements, 
gravity effects, significant persistent air turbulence, trans- 
lational coupling dynamics, and CMG desaturation using 
ducted fan based vector thrusting has been developed. 

An alternate control viewpoint is based on the linear 
dynamics in Figure 4. This perturbed dynamical system 
can be viewed as an linear parameter varying (LPV) sys- 
tem for control law design wherein the scheduled param- 
eters may include the trim parameters such as attitude 
and translational velocity. This viewpoint is particularly 
attractive because the control law analysis and synthe- 
sis is well understood and technologically viable (see for 
example [11], [12], [13], [14]). For further details on both 
control law approaches, see [8]. 

4 Flying Test Platform 

4.1 System overview 

The attitude control approach described above is tested 
using the vehicle pictured in Figure 8 as a CAD model. 
The vehicle is a platform levitated by four ducted fans 
arranged symmetrically about the 2 -body axis. Each fan 
duct is equipped with a parallel pair of thrust vectoring 
vanes, downstream from its fan. Each fan and pair of 
vanes is individually commanded by the control system. 
The 2 -body axis angular momentum bias is provided by a 
momentum wheel, powered separately from the levitation 
fans, and can be seen in the Figure at the center of the 
platform. Although the simplest and lightest momentum 
biasing mechanism is provided by augmenting the 2 -body 




axis moment of inertia of the propulsion system, momen- 
tum wheel is separately mechnized in order to facilitate 
tests in which different magnitudes of momentum bias 
are explored. For initial testing, the vehicle will operate 
without CMGs or reaction wheels. 

In its initial configuration, the fan motors and 
momentum wheel motor are Aveox F27-F5B, and JR 
DS8411 are the vane servos which provides 300 deg/s 
maximum deflection rate. Each of the four levitation 
fan motors draw two kilowatts at full thrust. Because 
of this, the vehicle will be fed power from an overhead 
tether whose tension is controlled to offset its weight. The 
tether will mount to the platform at the vehicle’s center 
of mass, to minimize tether-induced moment disturbances 
during flight. The dynamics of the tether are ignored in 
the next Section’s simulation model. The attitude sta- 
bilization software is hosted on a PC104 computer that 
is affixed to the platform, and which communicates with 
a ground-based dSpace real-time control system. This 
latter is used to generate trim and guidance commands. 
The platform hosts an IMU model MIDG II, which is 
supplemented by a machine vision subsystem that gives 
Earth-fixed position and attitude. More details of this 
system are given in [8]. 

Each fan produces roughly 10 lbf of thrust at full 
power, which is more than adequate to levitate the gross 
platform mass of approximately 29 lbm. The moments of 
inertia (in Kg-M 2 units) for this vehicle are I xx = 0.59, 
I yy = 0.58, I zz = 1.15, and /„■ « 0, Vi ^ j. 



Figure 8: The NASA Flying Test Platform. 


4.2 Bench Testing 

In order to generate control design models and perform 
detail dynamical simulations which are physically rele- 
vant, we first characterize by parameterizing the dynam- 
ics of the subsystems and total system from first princi- 


ples and then use a series of bench tests to determine these 
parameters from measured data. Specific necessary tasks 
to this end includes design and fabrication of specific test 
configurations, instrumentation, calibration, followed by 
model development and validation through estimation of 
parameters including noise levels and its associated fil- 
ters. 

Figure 9 shows the bench test configuration for an 
individual ducted fan and vane system. The focus in this 
bench test involves characterizing (i) Fan motor and vane 
dynamics, and (ii) Fan and vane aerodynamics. The for- 
mer basically involves the development of models to char- 
acterize fan speed and vane angle from their commands, 
while the latter modeling task can be anywhere from al- 
most trivial to near impossible if the unsteady aero effects 
are to be included with any level of confidence. Together 
they determine the actual net forces and torques gener- 
ated by control commands for fan speeds and vane deflec- 
tion angles. 


Figure 9: Bench test configurations and modeling. 

Initially we consider a static aero model for the con- 
trol forces and torques generated from fan speed and vane 
angle. For example, consider determining a steady aero 
model of net forces due to fan and vane. Figure 10 shows 
the measured net forces as a function of fan speed, for a 
fixed vane angle, while Figure 11 shows normalized mea- 
sured vertical forces as a function of van angle, at constant 
fan speeds. Figure 12 illustrates the measured unsteady 
effects in the vertical force component with time vary- 
ing fan speeds with fixed vane angle. The l-cr variation 
of 0.1 pound thrust in a single fan is not insignificant 
for this subscale platform, especially considering that we 
have not included additional factors such as vane effects 
and implicit bandwidth limitations in the fan motor. 

The above sample data were obtained experimentally, 
by testing the components in a ducted fan configuration, 
displayed in Figure 9. Static thrust vector force compo- 
nents were obtained from these tests, and are provided in 
[8] . Unsteady flow forces which are difficult to character- 
ize are not quantified, and do not appear in the simulation 
model described in the next Section. 
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Figure 10: Measured forces with fan speeds, at fixed 0 
degree vane angle. 


Figure 11: Measured and quadratic fit vertical force with 
vane angles, at different fan speeds. 


4.3 Simulation 

In this section we present simulation results based on full 
nonlinear models that are constructed from first princi- 
ples, bench test data refined model parameters, sampled 
and quantized measurements, and saturation limits on 
actuators. In particular, the dynamical models for the 
Ducted fan/vane, whose development are based on bench 
tests, includes asymmetrical effects which are difficult to 
visualize due to their 3-dimensional characteristics. This 
MATLAB- Simulink simulation model also incorporates 
models of servo dynamics that regulate fan speeds and 
vane angles. 

In all simulation cases which follows, the vanes and 
the fan speeds are actively used for stabilization and con- 
trol which reflects the current state-of-the art approach 
which is vector and differential thrusting, albeit, RC- 
quality indoor lab experiments. Initially, we focus our 
simulation on the effects of bias momentum in the stabi- 
lization and control performance of our flying test plat- 
form and do not include CMGs. We highlight these effects 
by comparing a system with significant level of bias mo- 
mentum (BM cases) with a corresponding system with 
no bias momentum (NMB). A difference in our simula- 
tion studies is to independently control all four fan speeds 
and their respective vane angles, as opposed to constrain- 
ing it to work in pairs and collectively, as in conventional 
helicopter control [18], [19] or practically all VTOL flying 
platforms [1]. 

4.3.1 Crosswind and turbulence model 

To simulate wind disturbances during hovering or low air- 
speed operations, crosswinds are assumed which impinges 
on a vertical cylinder attached to the upper part of the 
platform. The resulting aerodynamic drag on this cylin- 


drical column is used to simulate the crosswind and tur- 
bulence on the vehicle. The resulting disturbance model 
consists of forces and torques in the pitch and roll axis. 
Two sets of winds are simulated, namely, “light” winds, 
and “strong” winds. The “light” winds simulate “Calm 
to Gentle Breeze” with sample mean wind speeds of 7 
Knots and standard deviation of 5. However for “strong” 
winds, intended to simulate “Calm to Strong Breeze/Near 
Gale”, the simulated sample mean wind speed is 15 Knots 
with a standard deviation of 18, as seen in Figure 13. 

4.3.2 Controller Law Synthesis 

Based on LQG control theory, feedback controllers are 
designed to stabilize and hold the platform attitude and 
also track the translational velocity and angular velocity 
commands. In this simulated design study, feedback con- 
trollers for the BM and NBM cases are synthesized based 
on the linearized model about hovering trim, which is 
identical to the trim conditions at a given translational 
velocity if aerodynamics due to significant airspeeds are 
not included, as assumed in this study. Note that the 
trim conditions are necessarily different for both cases 
since aero-resistance fan drag torques in the case must be 
statically canceled by vane-deflected forces. 

The augmented linearized NFTP model with fan and 
servo models has 21 states consisting of 3 attitude pa- 
rameters, 3 linear momentum, 3 angular momentum, 8 
fan motor states, and 4 servo motor states. The 8 fan 
motor states includes the PI controllers which regulate 
each fan speed. The outputs of the augmented model 
are translational and angular velocities, both expressed 
in vehicle fixed coordinates. For a fair comparison be- 
tween BM and NBM cases, the same LQG weights were 
used in the controller design 
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Figure 13: Simulated strong wind. 
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Figure 12: Measured uncertainty in the vertical force with 

time varying fan speeds at a fixed vane angle. 4.3.3 Hovering 


Figure 14 shows the tradeoff between performance 
cost and control effort for a parameterized set of LQG 
controllers undergoing gust response during hover over a 
period of time Tf ina i. Specifically, the performance cost 
in terms of attitude and position errors are 




: 3 

E 

i=l:3 


[(ft-/3i) 2 + (&-&) 2 ]df (46) 


Open Loop Response to Light Winds. The open- 
loop attitude responses of the initially trimmed system for 
the two cases are shown in Figure 15. In the NBM case, 
the vehicle rapidly loses orientation and goes unstable, 
but in the BM case, the vehicle holds attitude inspite of 
a small gyric roll and pitch angle response to the light 
winds. However, the platform drifts off slowly without 
feedback compensation due to a steady wind component 
in the light winds. 


and the control effort is defined by a cumulative control 
effort (CCA) function as follows 

CCA(t) = [ (9l + 9l + 9l + 9l)dt (47) 
Jo 

The abscissa and ordinate in Figure 14 represents 
CCA(T fi na i ) and J(Tf ina i) respectively. In addition to 
quantifying the level of control activity, the above con- 
trol effort is intended to capture the level of control vane 
activity, which is seen as the primary source of control 
surface induced unsteady aerodynamics due to its vector 
thrusting role. Albeit limited, this tradeoff study indi- 
cates that controllers are significantly more effective for a 
system endowed with significant levels of bias momentum 
(asterisk line) in terms of actuator effort necessary for a 
given level of performance cost. The figure also shows the 
effects of a hypothetical model of unsteady aero depen- 
dence on actuator effort [8] , namely, the expectation that 
controlled performance can inadvertently degrade with 
overly active control surface activity as mostly seen in 
the NBM case (diamond line). However, in the following 
simulations, this hypothetical model of controller induced 
unsteady aerodynamics will not be included. 


Closed Loop Response to Light Winds. The 

closed-loop responses to light winds during hover for both 
cases are shown in Figure 16. Although both cases show 
stable response, the angular response for the BM case 
is clearly superior to the NBM case. Apparently, un- 
der these simulated light wind conditions, the angular re- 
sponse amplitudes for the NBM case appears to be small 
enough to be tolerated during an actual flight. The time 
history of the cumulative control effort (bottom right of 
Figure 16), shows that the case with BM actually requires 
less vane actuation than in the corresponding NBM case. 
As expected, the vane and fan actuator response did not 
saturate for both cases, under light winds. 

Closed Loop Response to Strong Winds. The 

closed-loop responses to strong winds during hover for 
both cases are shown in Figure 17. The controller for 
NBM case was judiciously scaled to mitigate actuator 
saturation. Both cases remains stable but the angular 
response for the BM case is clearly superior to the NBM 
case. However, the significantly larger angular response 
amplitudes for the NBM case will likely be unacceptable 
or even dangerous during an actual flight. The time his- 
tory of the cumulative control effort shows that the case 
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Figure 14: Performance cost (ordinate) vs actuator effort 
(abscissa) for an LQG set of controllers undergoing gust 
response during hover. 


Figure 15: Open loop response to light winds, 3-2-1 Euler 
angles. 


with BM actually requires less vane actuation although 
it gives far superior attitude hold than the corresponding 
NBM case. Some saturation in the vane and fan actuator 
resulted for the NBM case. 

To examine what happens if the controller for NBM 
case is not judiciously scaled to mitigate actuator sat- 
uration, resulting possibly to a more agressive feedback 
control, an additional simulation was done with the new 
control law. It turns out that the close loop system loses 
control after a few seconds. More specifically Figure 18 
shows the time histories of the vane angle and fan speed 
responses. The vane and fan actuators saturate resulting 
in loss of vehicle under strong winds. 

4.3.4 Command Tracking 

In this section, we outline a sample of simulation results of 
command tracking under strong winds, specifically trans- 
lational velocity commands while holding a trim attitude. 
It should be noted that the demonstrated inherent advan- 
tage of bias momentum in endowing directional stability 
to the platform do not necessarily improve its maneuver- 
ability. Hence, the simulation results in this subsection is 
intended to examine this issue quantitatively. 

For the following simulations, LQG based controllers 
are designed to track a trajectory path defined by a com- 
manded translational velocities about trim with zero an- 
gular velocity. The input commands are shaped to mit- 
igate the discontinuity in a step command. Figure 19 
shows the response to translational velocity tracking com- 
mands for both cases subject to strong winds. The accom- 
panying excursions in the angular velocities (and hence 
the vehicle attitude) are shown in Figure 20. It is seen 
that for both cases, the vehicle can successfully track 
translational velocity commands but BM case gives su- 


perior performance. 

The corresponding vane and fan actuator time histo- 
ries during tracking are shown in Figures 21 and 22, re- 
spectively. The case with BM clearly needs less vane and 
fan control activity but as observed earlier gives smaller 
velocity tracking errors. 

4.3.5 Robustness to payload variations 

Simulated Mass &: Inertia changes. To simulate 
variable payloads, we assume a step increase of 20plat- 
form. A spherical uniformly distributed mass is assumed 
dropped at (x=0.15,y=0,z=0) location on the platform at 
5 sec and is assumed to be dropped off from the platform 
at 20 sec. The inertia matrix is significantly changed as 
a result of the off-centered loading. This off-CG centered 
loading simulation is of particular interest due to its diffi- 
culty and its practical significance for a hovering platform 
where passengers or payloads can be loaded and unloaded 
with maximum flexibility without compromising vehicle 
stability while hovering. In this part of simulation, dis- 
turbances and measurement errors are not included to 
focus entirely on the effects of these payload changes. 

Closed loop response during hover. Figures 23 and 
24 show the position and attitude responses to mass and 
inertia changes for both cases. It is seen that the platform 
position and attitude excursions are significantly better 
in the BM case. 

Significant differences are also noted in their corre- 
sponding vane and fan actuator histories in Figure 25. It 
is seen that the fan speeds and vane angles oscillates sig- 
nificantly and rapidly in the NBM case resulting in the 
platform oscillations shown previously. Unfortunately, 
the actual consequences on flight performance due to this 
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Figure 18: Closed loop response to strong winds, more agressive control. 







Figure 16: Closed loop response to light winds. 


Figure 17: Closed loop response to strong winds. 


apparently more unsteady control vane activity for the 
NBM case, remains to be seen during an actual flight due 
to the lack of a reliable mathematical model for this class 
of unsteady aero. 

5 Conclusions 

A new stability and control technology, based on trying to 
mitigate control surface actuator induced unsteady aero 
effects in feedback control is proposed. Its ultimate goal is 
to enable statically unstable thrust levitated low-speed or 
hovering vehicles to operate more safely than is possible 
today, in turbulent conditions and under large payload 
variations. New results in trim, bias momentum sizing, 


and a control structure framework which integrates con- 
ventional vector thrusting with multiple wheel actuators 
have been developed. Simulation results indicate that 
generic VTOL platforms can be endowed with signifi- 
cant directional stability with a proper amount of bias 
momentum, similar to a dual-spin stabilized spacecraft 
but under a significantly different operating environment. 
The results also indicate that this open loop directional 
stability carries over to closed loop if the controllers are 
designed properly. Specifically, it was found that vehicle 
attitude robustness can be significantly improved (over 
a no-bias-momentum vehicle) during hovering, particu- 
larly under strong turbulent winds or significant payload 
variations. In addition, command tracking performance 
at low airspeeds can also be improved significantly par- 
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Figure 25: Actuator responses to step changes in payload. 
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Figure 19: Vehicle velocity response during trajectory Figure 20: Vehicle attitude response during trajectory 
tracking. tracking. 


ticularly under turbulent winds, despite the presence of 
somewhat complicated gyroscopic coupling. These per- 
formance improvements are due to control laws that are 
based on proper modeling of the system dynamics. These 
predicted improvements in robustness obtained from de- 
tailed simulation studies are expected to be more evident 
during actual flight test, with the physical presence of 
control induced unsteady aerodynamics. 
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